Bose Condensed Gas in Strong Disorder Potential With Arbitrary Correlation Length 
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We study the properties of a dilute Bose condensed gas at zero temperature in the presence of 
a strong random potential with arbitrary correlation length. Starting from the underlying Gross- 
Pitaevskii equation, we use the random phase approximation in order to get a closed integral equation 
for the averaged density distribution which allows the determination of the condensate and the 
superfiuid density. The obtained results generalize those of Huang and Meng (HM) to strong 
disorder. In particular, we find the critical value of the disorder strength, where the superfiuid 
phase disappears by a first-order phase transition. We show how this critical value changes as a 
function of the correlation length. 

PACS numbers: 03.75.Hh,03.75.Kk,05.30.-d 



I. INTRODUCTION 

An ultracold atom gas in the presence of disordered 
environments is becoming a subject of increasing experi- 
mental and theoretical research activities. Generally, one 
would like to understand how the condensation and the 
superfiuid properties of ultracold gases are influenced by 
a spatially random force on the atoms. In some experi- 
ments, the random potential is created by optical means 
to show its effects on the transport properties of a Bose 
gas P, H, In most others, however, one must rather 
face the reality of unavoidable external random forces 
which are induced either by the roughness of a dielectric 
surface Q , by the magnetic field along wires with current 
irregularities ||, or by different localized atomic species 
Furthermore, recent theoretical results on the im- 
pact of randomness on bosons in lattices are reviewed in 
Refs. 0,0. 

Most theoretical studies on a 3D disordered Bose gas 
are limited to calculations up to the second-order in the 
random potential. In these works, the weak disorder in- 
duces only small corrections to the condensate depletion, 
the superfiuid density [9(, the collective excitations and 
their damping (To| . and the condensation critical point 
fill ]. The extension to strong disorder has so far been 
analysed only numerically in Ref. [13] . More recently, 
an analytical mean- field study [HI], which takes into ac- 
count higher order corrections, has shown the possibility 
of having a transition from a superfiuid phase to a Bose 
glass phase where the spatial long-range correlations have 
completely disappeared. 

In this paper, we address the issue of the influence of 
strong disorder at zero temperature for a finite corre- 
lation length 0, [lj|. The random potential follows a 
Gaussian distribution and is said to be uncorrelated in 
the case where all Fourier components contribute equally 
to the randomness, while it is correlated when the influ- 
ence of the Fourier components falls off for wavenumbers 
larger the inverse correlation length £. As encoutered 
in experiments we choose a Lorentzian correlation 



function. Usually, this length appears to be much bigger 
than the healing length and thus affects the condensate 
properties. 

In order to simplify this physical problem, we assume 
that all particles occupy the same quantum state, for 
which the macroscopic wave-function obeys the Gross- 
Pitaevskii equation in the presence of an external spa- 
tially random force. In order to solve this stochastic non- 
linear differential equations, we apply the random phase 
approximation (RPA) [l6| and take the ensemble average 
over all possible realisations of the associated potential. 
In the clean case without a random potential, this gap- 
less and conserving approximation has been successfully 
used in the context of calculating of the collective excita- 
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tions at finite temperatures [171 . Il8j and in kinetic theory 
[l9| . In our case with disorder, we obtain the particle 
density distribution beyond the lowest order expansion 
in the random potential. With this we extend the sem- 
inal work of Ref. @ to strong disorder, which has the 
consequence that the superfiuid phase disappears by a 
first-order transition. We show how the critical value of 
the disorder strength, for which this transition occurs, 
depends on the correlation length. 



II. MOMENTUM DISTRIBUTION 

We start from the Gross-Pitaevskii equation at 
zero temperature written in Fourier space. Defin- 
ing the Fourier components according to iJ)(t) — 
Ek e ' k ' ra k/V / ^! we get in units with h = 1: 
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where the random potential C/ q follows a Gaussian distri- 
bution (UqUqi) = (5 q ,_ q 'i?(q)/V r . The quadratic ampli- 
tude is assumed to be Lorentzian i?(q) = i?/(l + £ 2 
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i.e. it is correlated below the wavenumber l/£. We 
assume that a macroscopic fraction of the condensate 
moves with a velocity k s /m. 

In order to derive a dynamic dielectric function for the 
total fluid, we define the bilinear combination 



a k «k+q • 



(2) 



It represents an excitation of momentum q created from 
a particle which releases its momentum from k + q to 
k. In particular, it allows the definition of the Fourier 
component of the density fluctuation: p q = pk,q = 
p*_q. The dynamic evolution of ||3J| following from JT]) is 
given by: 



d 

'at"*- 



q' 



= ( e k+q — £k)Pk,q + 

+ |^X! Q: k' Q; k'+q')( Q: kak+q-q' 



(3) 

"k+q'^k+q) • 



The technical details for implementing the random phase 
approximation in the clean case can be found either in 
Ref. [16| or in Ref. [19| where a finite-temperature Bose 
gas has been considered with pk,q being a quantum op- 
erator. Here, we shall briefly repeat this procedure by 
considering random variables instead of operators. 

We treat the quartic terms in ([3]) within a factorisa- 
tion procedure and use the property that (pk.q) = for 
q 7^ in an homogeneous gas. For any quartic term, we 
approximate 

"ki^ka^ksaki - (a kl a k2 ak 3 a!k4) — 

+ (|Okx| 2 )^i,k3 a ka a ki 

+ (l a k 2 | 2 )<5k 2 ,k 4 (l + £ki,k 2 £k 3 ,k 4 )a k:L a!k3 

+ (|ak 1 | 2 )^ki,k 4 (l - <5k!,k 2 - <5k3,k 4 )ak 2 a k3 

+ (|ak 2 | 2 )<5k 2 ,k3(l - ^ki,k 2 - <5 k3 ,k 4 )a kl a k4 (4) 

which avoids double counting. For example, in case of 
ki = k2 and k3 ^= k4 the approximation reduces to two 
terms only which is important for contributions involving 
the macroscopic component k s of the condensate. Since 
the average over the quartic term in Eq. (j4|) applied in (|3|) 
involves components with the total transfer momentum 
q 7^ 0, it will not contribute for an homogeneous gas. 

Through this procedure, the RPA keeps among all 
terms those combinations involving products of off- 
diagonal terms p k ',q and averaged diagonal ones n k " = 
(l a k"| 2 ) f° r an possible values of k' and k", and neglects 
all others combinations. Therefore, we remove contribu- 
tions which are bilinear in p k ' !q ' for q' ^ q, 0. In this 
way, we obtain the linear integral equation for q ^ 0: 
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Here n' k = (1 — <5 ks)k )n k refers to the disordered part 
of the condensate which consists of all its part which are 
not at the wavenumber k s . Note at this stage that p kjq is 
still a random variable. Thus, we should still take the en- 
semble average over any non-vanishing combination like 
(J7_ q p k! q) and (pk' ,-qPk,q) and solve the resulting equa- 
tions. Equivalently, here we directly solve Eq. ((5]) for p k q 
and perform the disorder average at a later stage. 

In order to make the link with the dielectric for- 
malism, let us assume for the moment that the poten- 
tial has a temporal dependence of the form Uq(t) = 
cxp(— iojt)Uq jW . Then the solution is of the form 
Pk,q(i) = exp(— MJt)pk,q,w as well. Under these condi- 
tions, we find a solution: 
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which is similar to Ref. [19]. The Fourier component of 
the density fluctuations can therefore be written in the 
form: 



with the susceptibility 
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The dynamic dielectric function /C(q, to) defined by 
for the total fluid can be decomposed into [l9[ 



/C(q,w) = /C„(q,w) 



2gn ks q 2 
V m 
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where we obtain for the disordered part of the fluid: 



/c„(q,^) = i- vE 



'k+q 



v — u . 



' m 2m 



(10) 



(cf. also [ijJUO])- In this way, we recover the same ex- 
pression as in Ref. [13, H3] for the susceptibility function 

x(q,w). 

In the special case of a time-independent potential 
Uq,ui — Uq6 U Q which is of interest here, the solution © 
can be used to define a self-consistent relation that allows 
the calculation of n k . At first, combining Eqs. (fB|)— (fTUj) . 
we get the response function for a particle of the conden- 
sate to be excited with momentum q: 
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For cj = 0, the screening factor AC(q) denned by (fTTj) for 
any external force acting on the condensate particles is: 



K(q) 



V m ' T V 2m / 
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In the simple special case n k = 0, i.e. /C„(q, 0) = 1 and 
k s = 0, we obtain more simply: 



£(q) = 



1 



1 + (4m 5 n /T/)/q 2 



(13) 



This formula shows that the random force acting on the 
condensate particles is screened for momenta below the 
inverse of the healing length which plays here the role of 
the Debye length for the condensate. 

Let us now notice that: (p£ s q Pk s ,q) = 
(|akj 2 |ak s +q| 2 ) = «k B " ks+q for q ^ whereas the 
fluctuations of the macroscopic k s component around 
the average are small in the thermodynamic limit. Using 
this relation in (fTTj) and (| 1 2 [) . we arrive at the following 
non-linear integral equation for n'^^ : 
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where = ^/c^q 2 + (q 2 /2™) 2 denotes the Bogoliubov 

excitation energy and cb = y gn^/mV is the sound ve- 
locity in the absence of disorder. The nonlinearity comes 
from the fact that the dieletric function /C„(q, 0) depends 
via Eq. (flT))) on n ks+q - Let us first consider the case of 
weak disorder where we can approximate /C„(q, 0) ~ 1 
and recover the second-order result for the disordered 
components of the condensate [§] : 



^(q)(£k s 



(15) 
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FIG. 1: Clean part of the condensate fraction n as a func- 
tion of R* in the case of uncorrelated disorder in the RPA 
model (full curve) and in the HM model (dotted curve) and 
the corresponding superfluid fraction n s (dot-dashed curve in 
RPA and dashed curve in HM model). 




FIG. 2: Same as for Fig.l but for correlated disorder with 

r = i. 



This formula is regular provided the Landau stability cri- In the limit of sma11 velocit y k */ m > we § et the expression: 



terione q > |k s .q/m| is satisfied. An increase of k s would 
increase the disordered part of the fluid until a singular- 
ity is reached leading to an instability. For k s = 0, we 
recover the weak-disorder result for this part 
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Noticing that /C„(q, 0) with the solution n' k is an even 
function of k s , we deduce from Eq. (fT4|) : 



The superfluid fraction n s is dehned as the part of the 
fluid which is moving at a velocity k s /m and is found by 
calculating the total momentum of the gas: 



n - n s = - (n - n ) . 



(19) 



P = y^kn k = k s n s 



(17) 



This relation is identical to the one in Ref. [PJ but is here 
found to remain valid in the RPA for the more general 
case of strong disorder. 
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FIG. 3: Same as for Fig.l but for correlated disorder with 

e = 10. 
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FIG. 4: Critical value for the disorder intensity R* as a func- 
tion of the correlation length f*. 

III. PHASE TRANSITION 

Equation (fT4"| is solved numerically by an iterative pro- 
cedure as a function of the reduced dimensionless param- 
eters for the disorder strength R* — Rm?/ 2 / ^An 2 gn/V 

and the correlation length £* = y^Amgn/V^. We have 
looked for the stable solution that minimizes the total 
energy of the system. The homogeneous part of the 



condensate no and the superfluid fraction n s are plot- 
ted in Figs. 1-3 as a function of the reduced disorder 
strength R* and are compared with the results obtained 
by HM . The transition is determined as the highest value 
R* = R* for which the condensate fraction uq is non- 
zero so that spatial coherence is preserved over the entire 
space. This point corresponds to the situation where no 
has an infinite derivative with respect to R* . As a conse- 
quence of relation (fl~9f , also the superfluid density n s has 
there an infinite derivative. For uncorrelated disorder i.e. 
£* = 0, we notice only a tiny difference with respect to 
the second-order theory, until we reach the critical value 
when the disordered part of the condensate is about 30 
percent. For higher values of £* this difference becomes 
more pronounced, but a smaller fraction of the disordered 
part of the fluid is needed in order to achieve the transi- 
tion. Fig. 4 shows that the critical value R* c increases as a 
function of the coherence length which is understandable 
from the fact that for larger £* the high spatial frequency 
components of the disordered part of the condensate are 
smaller (cf. Eq. (fT4")) ). 



IV. CONCLUSIONS 

The random phase approximation has been applied to 
the Gross-Pitaevskii equation in the presence of a ran- 
dom potential in order to describe a strongly disordered 
Bose gas at zero temperature. This approximation goes 
beyond a previous second-order calculation and predicts 
a first-order phase transition from a superfluid phase to 
a non-superfluid phase. The critical value of the dis- 
order intensity for this transition depends strongly on 
the correlation length. Nevertheless, our model fails to 
describe the properties of the non-superfluid phase. A 
possible explanation is that the assumption of a unique 
wave function for any particle excludes the possibility of 
having fragmented condensates for strong disorder that 
could be necessary in a such a phase. 
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